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Analytical expressions are obtained for the bulk density of monodispersed powder systems, taking account of 
the properties of the particles. The structural and percolation approaches as well as taking account of the 
boundary layer thickness substantially simplify the evaluation of the properties of powder filled compositions. 
It is shown that the computed data are good agreement with the experimental data. 
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Different structural and percolation models as well as 
different fractal methods are used to describe powder materi¬ 
als which are nonuniform on macroscopic and microscopic 
scales, specifically, ceramics [1 - 7]. 

To evaluate the properties of porous and composite 
mono- and bi-dispersed materials it is necessary to take ac¬ 
count of the decrease of the bulk density of the filler V r as the 
diameter of the particles decreases or the specific surface 
area increases, since it changes the characteristic of the sys¬ 
tem. For porous materials V r = 1 -P r (P r is the bulk poro¬ 
sity). We note that the standard volume concentration of free 
smooth spherical particles in bulk or their random packing 
density V r0 = 1 - P r() is 0.637. 

Let the clusters (compact accumulations of particles) be 
spherical on average and let the order be different depending 
on the porosity of the system and in accordance with the 
fractal approach, i.e., the clusters can be the particles of a 
cluster system of a higher order. Likewise, let the small parti¬ 
cles be distributed in the gaps between larger particles (or 
clusters of these particles) with random packing density V r . 
In systems with elongated particles the relative bulk density 
V r can be expressed as V r = 1.5 + (L/D) 9 ^/ 11, where D and 
L are, respectively, the diameter and length of the elongated 
particles [5]. 

In [5] a relation that gives a good fit to experimental data 
was proposed for bulk density of monodispersed systems: 
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Taking account of the shape and roughness of the parti¬ 
cles according to the specific surface area of the particles 
D = 6/S 0 the expression (1) can be written in the form 

F r =l-0.1(6/5 0 )- 019 
for 6000 m“ 1 < S 0 < 9.23 x 10 5 m- 1 ; 

V r = 0.63 for S 0 < 6000 m “ *; 

V r = 0.97 for S Q > 9.23 x 10 5 m-‘. 

However, it is useful to obtain analytic relations for the 
bulk density of monodispersed system taking account of the 
properties of the particles. 

For this we shall determine the content of the small and 
large particles in a bi-dispersed system with maximum den¬ 
sity, i.e., as D/D 0 —» 0 the maximum relative packing density 
of the particles in the system F max = 0.637 + 0.363 x 0.637 = 
0.868. Then the large-particle fraction in unit volume of 
the material with maximum particle size is C sizemax = 
0.637/0.868 = 0.734. 

Let us evaluate the minimum ratio of the sizes (diame¬ 
ters) of the large and small articles or clusters under the con¬ 
dition that the small particles (clusters) fall between the large 
particles. In [4] it is established that for ordinary dispersed 
systems it is sufficient to consider three-layer clusters; then 
the distance between the center of a cluster to the center of 
the particles of the outer layer is /= 1.82D (i.e. the average 
diameter Z), + { of the cluster formed is 2/ = 3.64 D t ). 

We assume that the degree of clusterization of the system 
increases with decreasing particle size (because of an in¬ 
crease of the interparticle interaction forces). According to 
the experimental data, for particles larger than 1 mm the de¬ 
gree of clusterization is zero and for dust particles (size close 
to zero) it can be regarded as maximum and equal to 0.734. 
This corresponds to the volume ratio of large and small 
phases in the highest density packing, studied above, of a 
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Fig. 1. Model of a two-dimensional cluster of first (a) and se¬ 
cond ( b ) orders. 

bi-dispersed system. Figure 1 shows a model of a two-di¬ 
mensional (for clarity) cluster. It is evident that it is almost of 
a fractal character (i.e., as the scale increases, like but 
smaller elements of a similar structure appear). 

Then the following expression can be written for the 
clusterization A of a monodispersed system of particles of 
clusters: 

A = 0.734(1 - D/D 0 ) or A = 0.734(1 - 6000/S 0 ). (2) 

By interpolation from the computational data of [5] we 
obtain for random monodispersed packing (average coordi¬ 
nation number with bulk density Z r = 7.3) a ratio of the 
intercluster to cluster porosity equal to 0.93. For this reason 
it can be assumed that with cluster formation the density of 
the clusters increases by 1/0.93 and the intercluster density 
does not change. 

It is convenient to start the analysis with the smallest 
first-order clusters and then analyze larger clusters — se¬ 
cond-order and so on. In a system containing the smallest 
clusters of size D, we take account essentially only of its 
nonuniformity (ratio of the intercluster-to-cluster porosity). 
We have for the density of such a system 

V 1 = 0.637[(1 -A)+A/ 0.93]. (3) 

For systems containing clusters of size D x and D 2 = 
3.64Z) 1 = 3.64 2 Z), we take account of the intracluster poro¬ 
sity. The result is 

V 2= V 1 [(1-A) + AV i/0.93]. 

We obtain a similar result for a system containing clus¬ 
ters of size D u D 2 , D { (D n _ i = 0.001/3.64*): 

V i =V i _ l [(\-A) + AV i _ i/0.93]; (4) 

V r =V n =V n _ x [{\-A)+AV n _j0.9^ 
for D 0 = 10“ 3 m >D > 1.56x 10 _6 m; 

u 5 I 

V r = 0.63 for D > 10 “ 3 m; 

V. = 0.97 for D < 1.56x 10 -6 m. 

' 9 


Pr 



Fig. 2. Bulk density D r of dispersed systems versus the average size 
D of the constituent particles: 1 ) fit of the function (1); 2) experi¬ 
mental data (O); 3 ) calculations using the expressions (2) - (5). 


The reference values of the cluster sizes can be found 
from the relation D n _ f = 0.001/(3.64*): D n = Z) 0 = 0.001 mm, 
D n _i = 0.000275 mm, D n _ 2 = 0.0000755 mm, D n _ 3 = 
0.0000207 mm, D n _ 4 = 0.0000057 mm. The lower limit of 
the particle size is D mm = 1.56 x 10“ 6 mm, since calcula¬ 
tions of V r of the system give 0.0298, which corresponds to 
the maximum experimental values of the porosity of powder 
systems. 

Computational Algorithm. We shall determine the re¬ 
quired number of computational steps n (2 <n>6) with 
the aid of the auxiliary value n mx = 1 + (1 -D/D 0 )/ 3640, 
rounding the value obtained to the next highest integer. 

We calculate V r with the required number of computa¬ 
tional steps using the relations (2) - (5). 

Figure 2 displays the bulk porosity of powder systems as 
calculated from the expressions (2) - (5). The figure also 
shows the experimental data and the fit of function (1). Evi¬ 
dently, the calculations based on the expressions (2) - (5) are 
in good agreement with the experimental data and the func¬ 
tion (1). In contrast to the latter the proposed method reflects 
the particle sizes and the degree of clusterization as well as 
the fact that the intracluster porosity in the system is less than 
the intercluster porosity. 

The large variance of the experimental data is explained 
by the use of particles with different shape and roughness. To 
refine the calculations the particle size can be taken as the 
value calculated according to the average specific surface 
area of the particles D = 6/ S 0 . 

It is of interest to evaluate the fractal size of the cluster 
model described above. For a random system of spheres with 
random packing V r = 0.637 and D n = 1 we obtain 

^ = 6F„/(3.l4Dj); 

where N n is the number of clusters per unit volume of the 
materials with n computational steps; V n is the relative pack¬ 
ing density of particles in a system with n computational 
steps, D n . = D /3.64 = 0.2747T) , with density 0.637 in the 
space between the spheres D n and so forth. If one sphere 
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Fig. 3. Iinternal coefficient of friction jlx of powders as a function of 
the average diameter of the particles according the experimental 
data of [3]: O) A1 2 0 3 ; •) Si 3 N 4 ; O) Zr0 2 ; A) MgO; □) Al(OH) 3 ; 
V) AIN; +) TiN; solid curve) according to the expressions (6). 


Fig. 4. Modulus of elasticity of filled polymer systems: 1 ) fraction 
120-400 pm with bulk density 59.8% (quartz); 2) fraction 
5-10 pm with bulk density 49.0% (quartz); 3) fraction 0.001 - 
0.004 pm with bulk density 10.1% (aerosil-175); 1-3) calculation 
using the expressions (8) and (9); points) experimental data of [8]. 


(N n = 1) is used as the initial sphere, then according to the 
theory of fractals [6] we have 

D, = D n Z' = Z", 

where Z is the order of the self-similar relation. 

According to the Hausdorff- Besikovitch relation taking 
account of the fact that Z= 1/3.64 and the number of parti¬ 
cles (or lower-order clusters) in a cluster p = 28 [6], we ob¬ 
tain for the dimension d of a disperse system 

d= In p /In (1 /Z) = In 28/ln 3.64 = 2.58. 

Therefore, the fractal character of the powder system is 
confirmed. 

For systems with average particle size less than 6.5 pm, 
where there is no correlation between the porosity and the 
average particle size in the system, it is convenient to use the 
internal coefficient of friction p 


When using instead of the experimental data on V the ef¬ 
fective values of the concentration V eff according to the ex¬ 
pression (7) the experimentally determined threshold (criti¬ 
cal) values V CY will not contradict percolation theory, and 
then other structural characteristics and technological partic¬ 
ulars of obtaining a polymer filled composite material can be 
analyzed. 

We shall examine the percolation estimate of the modu¬ 
lus of elasticity of polymer systems with a filler based on sil¬ 
icon oxide (granite, quartz, aerosil), making it possible to 
take account of the effect of the filler density in a free-flll 
state V r on its critical density V cr [5]: 


E = (E 0 



^ V(D + cl) i / D 3 - V cr ' 
yV r (D + h) 3 / D 3 - V C[ j 



H =0.5 + 0.03(£>-6.5) for 1(T 6 <£><6.5x 10~ 6 m; 1 
\x = 14 + (1 -D) 4 exp (1—£>)+ 0.5165 for D < l(T 6 m.j 

Figure 3 shows the dependence of the internal friction 
coefficient on the average particle diameter according to the 
expressions (6) and the experimental data of [3]. It is evident 
that the computed and experimental data are in good agree¬ 
ment with one another. 

When evaluating the properties of composite materials, 
aside from the decrease of the relative bulk density of the fill 
V r with decreasing particle diameter, it is also necessary to 
take account of the higher mechanical properties of a poly¬ 
mer (boundary) layer with thickness hj of the order of 
0.15 pm adjoining the surface of the filler particles or the in¬ 
crease of the effective particle size by this thickness. As a re¬ 
sult the effective volume concentration of the fill V eff can be 
expressed as 

V eS =V r (D + 2h f f/D\ (!) 

where V is the relative density of the fill and D is the average 
(over the squared radius) size of the fill particles. 


where E , E 0 , and E are, respectively, the modulus of elasti¬ 
city of the filled polymer, the filler powder with the bulk 
density, and the polymer; V cr is the threshold density of 
monodispersed powder or fiber systems, which is related 
with the bulk density V r by the following expression [5]: 

V cr = V r ls /2.82 +0.003. (9) 

The expressions (8) and (9) with E 0 ^ = 1.1 x 10 4 
(2.2 x 10 3 ) MPa determined from the data of [8] were used 
to calculate the modulus of elasticity of different fractions, 
filled with mineral components (quartz, granite, aerosol), 
with different density of the polymer systems (based on 
furfural acetone resin). The computational results are pre¬ 
sented in Fig. 4. As one can see, the calculations are in good 
agreement with the data. 

In summary, the structural and percolation approaches 
make it possible to simplify considerably the evaluation of 
the properties of powder systems and filled polymer compo¬ 
sites. The packing ratios of the particles of the mineral fill 
and the thickness of the boundary layers must be taken into 
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account. The expressions (7) - (9) can be used, for example, 
to optimize the compositions and properties of composites 
containing an abrasive filler and powder binder consisting of 
glass in the fabrication of the ceramic abrasive articles. 

The approach proposed here makes it possible to refine 
the modeling of the properties of porous powder materials. 
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